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Abstract. 

^ c"| In a Pulsar Wind Nebula (PWN), the lifetime of inverse Compton (IC) emitting e~ exceeds the 

Q_[ lifetime of its progenitor pulsar (as well as its shell-type remnant), but it also exceeds the age of those 

that emit via synchrotron radiation. Therefore, during its evolution, the PWN can remain bright in 
IC so that its GeV-TeV gamma-ray flux remains high for timescales much larger (for 10 5 — 10 6 
yrs) than the pulsar lifetime and the X-ray PWN lifetime. In this scenario, the magnetic field in 
the cavity induced by the wind of the progenitor star plays a crucial role. This scenario is in line 
with the discovery of several unidentified or "dark" sources in the TeV gamma-ray band without 
X-ray counterparts; and it is also finding confirmation in the recent discoveries at GeV gamma rays. 
Moreover, these consequences could be also important for reinterpreting the detection of starburst 
galaxies in the TeV gamma-ray band when considering a leptonic origin of the gamma-ray signal. 
Both theoretical aspects and their observational proofs will be discussed, as well as the first results 
£^ . of our new modeling code. 
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ANCIENT PULSAR WIND NEBULAE 



^ ■ Although they are often detected as non-thermal X-ray sources, evolved PWNe can indeed 

lead to a bright y-ray source without any counterpart [1] [2] [3]. The key issue is that the 
low-energy synchrotron emission, depends on the internal PWN magnetic field [4] which 
may vary as a function of time, following Xe <*■ t 2a if B(t) t~ a , where %e is the syn- 
chrotron emitting lifetime of TeV y-ray emitting leptons with energy E and a is the power- 
law index of the decay of the average nebular field strength, whereas the VHE emission 
depends on the background radiation fields, which are constant on time-scales relevant for 
PWN evolution. MHD simulations of composite SNRs [5] find that a — 1.3 until the pas- 
sage of the reverse shock (RS); after the field decay should continue for much longer time 
since expansion continues. From hydrodynamic simulations, an expression for the time 

of the RS passage has been given [7] as T R = 1 Okyr ( w J^ cm ^ )~H ^fc ) * ( ^ ) ~ I 

where the I st term is the density of the ISM, the 2 nd the ejecta mass during the SN explo- 
sion and the 3 rd the SN blast wave energy. The stellar wind of a high- mass star can blow 



a cavity around the progenitor star [6] with relatively low ISM density, so that Tr>> 10 
kyr. In such a case B(t) can decay as t~ a , until the field is low enough for the X-ray flux 
to drop below the typical sensitivity levels. As a result, in a scenario where the magnetic 
field decays as a function of time, the synchrotron emission will also fade as the PWN 
evolves. The reduced synchrotron losses for high-energy e~ for such a scenario will then 
lead to increased lifetimes for these leptonic particles. For timescales shorter than the 
inverse-Compton lifetime of the e~ ( tjc 1.2 x 10 6 (£ e /lTeV) _1 years), this will result 
in an accumulation of VHE e~ (e.g. seen in HESS J1825-137) which will also lead to an 
increased y-ray production via IC. To summarize, during their evolution PWN may appear 
as y-ray sources with only very faint low-energy counterparts and this may represent a vi- 
able model for many unidentified TeV sources. This effect can be clearly modeled, tested 
on young PWNe such as the Crab or G2 1.5-0. 9, and finally applied to some candidate 
ancient PWNe such as HESS J1427-608 and HESS J1507-622. 

THE NEW CODE 

In this code (work in progress) we adopted the time evolution of the lepton spectrum of 
[8] (eq. 7 and 11), and the source function for the system Q(E e ,t) given by [9] 

n(F ri = J Qr(t)(E e /E b y Pl if E mia <E e <E b 
^ " ' I Gx(0 {E e /E b )-» if E b <E e < E max 

where the normalisation factors Qi(t) are determined using the prescription that the 
energy in the lepton spectrum should be some fraction of the spin-down luminosity L(t): 
J Qi(E e ,t)E e dE e = r}iL(t). Note that the two normalisation factors are solved indepen- 
dently, with the only requirement that r} t + rf x < 1 . The time evolution of the spin- 
down luminosity is given by L(t) — j^fp^i, where T is the characteristic spin-down 
time-scale of the pulsar. The synchrotron T syn and escape T esc time-scales can respec- 

tively be approximated by [8] T syn (V) « 1.25 x 10 7 ( l -^y ) ( 1 gf^ ) Y r an d *esc(t) f» 

Since a = 1.3 until the passage of the RS, the evolution of the radius of the PWN can 
be approximated by 



r (A -I r o( ? Ao) L4 if ? < ?rs 
^"l r (t/t ) - 2 if t>t IS 



where ro is a normalisation constant and ? rs = 5000 yr the time needed for the RS of the 
supernova shell remnant to reach the PWN. For the PWN radius we assume a power-law 
expansion with index 1 .4 before the passage of the RS and index 0.2 after it. 



Testing the new code on young PWNe (G21.5-0.9) 



In order to test the consistency of our code we tested it on the Crab nebula and on 
G21. 5-0.9, obtaining compatible results with [8] and with its previous version [1] [3]. The 
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FIGURE 1. The SED of G2 1.5-0. 9 shows that our new model can describe well the spectral features of 
young PWNe. The data are taken from [3]. 
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FIGURE 2. left: SED of HESS J1507-622 modeled with the new code; the radio information is taken 
from [11], the ferm/-LAT points from [12], the Chandra and H.E.S.S. points from [2]. right. SED of 
HESS J1427-608 modeled with the new code; data are taken from [13]. 



derived value for the magnetic field of G2 1.5-0. 9 is similar to the value estimated for the 
950 yr old Crab Nebula. Other parameters derived are listed in Tab. 1. The results, shown 
in Fig. 1 and summarized in Tab. 1, demonstrate that the new code is describing well the 
spectral energy distribution (SED) of G21. 5-0.9. 

In order to fit the radio and X-ray synchrotron emission, it is necessary for the e 
spectrum to have a discontinuity between the two parts of the broken power-law. A 
similar feature has been discussed by [10] to account for the emission from the Vela 
PWN. Despite the different parent e~ spectrum, the results are compatible with [3]. 

The code describes well also Unidentified sources 

After testing the validity of the code on young, well-known PWNe, we test it on 
Unidentified high-energy y-ray sources and especially on the so called "dark sources", 
i.e. on the y-ray sources that, even after deep multi-wavelength observations, lack any 
plausible X-ray counterparts (such as HESS J1507-622, HESS J1427-608, HESS J1708- 
410 and HESS J1616-508). Here we show the results on HESS J1507-622 and HESS 
J1427-608. The results are summarized in Fig. 2 and in Tab. 1, showing that our new code 
describes the SEDs also of these 2 "dark sources" very well, strengthening the idea that 
we are dealing with ancient PWNe. 

It is important to underline that for HESS J 1507-622 (the best example for this kind of 



TABLE 1. Model parameters for the three sources discussed in the text. 



G21.5-0.9 HESS J1507-622 HESS J1427-608 



Lo (10 J7 erg/s) 

r age (yr) 

Distance (pc) 
T(yr) 

z?(r age ) (uG) 
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4800 

250 
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0.06 
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0.05 
6 x 10 2 



£max (TeV) 
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E b (TeV) 



4 



5 x 10~ 2 
0.6 



sources [2] and the most challenging for the model, especially given the very low radio 
upper limit) the age (15 kys) is perfectly compatible with the estimation derived in [3], 
as is the distance (6 kpc) with the lower limit of [2]. Moreover, since there is no pulsar 
detected in the region of HESS J1507-622 or HESS J1427-608 (as expected for such old 
PWNe that the pulsar could have already spun down below the sensitivity of the current 
instruments [5]), an arbitrary choice of the initial conditions for the pulsar parameters 
is required; moreover, in order to fit the data for HESS J 1507-622, we need a current 
magnetic field of about 0.3 /lG. This value is compatible with the upper limit of 0.5 
/iG derived in [2], however we should underline that this value is smaller than the ISM 
magnetic field (that would raise again the question of section 3 of [3]). 



Our new code represents an useful tool to describe PWNe, since it could describe so far 
both their young and aged phases. The preliminary results on the three sources shown 
here are summarized in Tab. 1. Together with the new high-energy observations [12], it 
is strengthening the idea that VHE unidentified sources can indeed explained as ancient 
PWNe. Regarding the used model, we notice that: (1) the e~ energy spectrum needs to 
be described by a discontinuous broken power-law in order to simultaneously fit all the 
MWL data; and (2) the fraction of spin-down power needed for the higher-energy part 
of the e~ spectrum is found to be ~10 times smaller than the energy needed for the 
low-energy part of the spectrum. The consequences on Galactic and extragalactic (e.g. 
Starburst galaxies [14]) astrophysics, and consequently on CRs studies, are remarkable 
and already discussed in [2] [3] [5]. 



SUMMARY AND CONCLUSIONS 
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